Hole spin relaxation in p-type GaAs quantum wires investigated by numerically 
solving fully microscopic kinetic spin Bloch equations 
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We investigate the spin relaxation of p-type GaAs quantum wires by numerically solving the fully 
microscopic kinetic spin Bloch equations. We find that the quantum-wire size influences the spin 
relaxation time effectively by modulating the energy spectrum and the heavy-hole-light-hole mixing 
of wire states. The effects of quantum-wire size, temperature, hole density, and initial polarization 
are investigated in detail. We show that, depending on the situation, the spin relaxation time can 
either increase or decrease with hole density. Due to the different subband structure and effects 
arising from spin-orbit coupling, many spin-relaxation properties are quite different from those of 
holes in the bulk or in quantum wells, and the inter-subband scattering makes a marked contribution 
to the spin relaxation. 

PACS numbers: 72.25.Rb, 73.21.Hb, 71.10.-w 



I. INTRODUCTION 

Spintronics is continuing to attract considerable in- 
terest because of its potential application to informa- 
tion technology.— Several spintronics devices have been 
proposed that manipulate spin via spin-orbit coupling 
(SOC). 2,3,4 In recent years, progress in nanofabrication 
and growth techniques has made it possible to produce 
high-quality quantum wires (QWRs) and investigate spin 
physics in semiconductor nanostructures^&i The energy 
spectrum of QWR systems with strong SOC has been 
extensively studied ^Ol^d 3 ^^^^^* 2 ! It is well 
known that, even in the absence of an external magnetic 
field, the RashbaS and Dresselhaus^ SOCs lift the spin 
degeneracy in wire subbands at nonzero wave vectors. 
The subband structure for quantum-confined valence- 
band states is particularly interesting since they are sub- 
jcct to an especially strong §00^^17,18,19,20,21,22 

As a long spin relaxation time (SRT) is desir- 
able for the operation of spintronic devices, many 
investigations have been performed to better un- 
derstand the electron spin relaxation in quantum 
ctrtrcturec , 23,24,25,26,27,28,29,30,31,32,33,34,35,36,37, 38,39,40,41,42 
e.g., using the single-particle mode l 23 ' 25 i 26 ' 31 i 32 ' 33 i 34 or 
Monte-Carlo simulations i 29 i 30 ' 31 However, it was shown 
by Wu et alS^^^^^- that the single-particle approach 
is inadequate in accounting for the spin relaxation. A 
fully microscopic kinetic spin Bloch equation (KSBE) 
theory, which takes full account of the inhomogeneous 
broadening from the Dresselhaus and/or Rashba SOC 
and the effect of scattering, has been developed to study 
spin relaxation i 35 ' 36 ^! 3 ^ 9 .^ . Cheng et al. applied this 
approach (excluding the Coulomb scattering) to study 
electron spin relaxation in QWR systems and showed the 
feasibility of manipulating spin decoherence4£ Investiga- 
tions of spin relaxation of holes in QWRs are relatively 



limited^ even though knowledge of hole spin relaxation 
in p-type QWRs is important for assessing the feasibility 
of hole-based spintronic devices^ 7 - The spin-relaxation 
mechanism in hole QWRs can be expected to be quite 
different from that in electron systems, and 2D or bulk 
hole systems, due to the strong SOC and the complex 
wire-subband structure. These effects have not been 
addressed previously. 

In this paper, we investigate hole spin relaxation in 
a p-doped (001) GaAs QWR. An idealized system of 
quantum wire with rectangular confinements and hard 
wall potential is considered in our calculation. First, we 
obtain the subband structure by diagonalizing the hole 
Hamiltonian including the quantum confinement. Here 
the light-hole (LH) admixture is dominant in the low- 
est spin-split subband, but the heavy-hole (HH) admix- 
ture becomes also important in higher subbands due to 
the strong HH-LH mixing. Then we investigate the time 
evolution of holes by numerically solving the fully micro- 
scopic KSBEs in the obtained subbands, with all the scat- 
tering, particularly the Coulomb scattering, explicitly in- 
cluded. We find that the QWR size influences the SRT 
^fi^iMvcly because the SOC and the subband structure in 
QWRs depend strongly on the confinement. When the 
QWR size increases, the lowest spin-split subband and 
the second-lowest spin-split subband will get very close 
to each other at an anticrossing point. If the anticrossing 
is close to the Fermi surface, the contribution from spin- 
flip scattering reaches a maximum and, correspondingly 
the SRT will reach a minimum. Moreover, we show that 
the dependence of the SRT on confinement size in QWRs 
behaves oppositely to the trend found in quantum wells. 
It is also found that, when the QWR size is very small, 
the SRT can either increase or decrease with hole density, 
depending on the spin mixing of the subbands. However, 
the behavior of holes in QWRs where the SRT increases 
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or decreases with hole density is quite different from the 
one of LHs in quantum wells with small well width. 45 
These features originate from the subband structure of 
the QWRs and the spin mixing which give rise to the 
spin-flip scattering. The spin mixing and inter-subband 
scattering are modulated more dramatically in QWRs by 
changing the hole distribution in different subbands. We 
also investigate the effects of temperature and initial spin 
polarization, showing that the inter-subband scattering 
and the Coulomb Hartree-Fock contribution can make a 
marked contribution to the spin relaxation. 

This paper is organized as follows: In Sec. II we set 
up our model and the KSBEs. Our numerical results are 
presented in Sec. III. We conclude in Sec. IV. 



II. MODEL AND KSBE 

Our investigation considers a rectangular p-doped 
(001) GaAs QWR confined in both x and y directions 
as shown schematically in Fig. 1. The potential height 
of the barrier layer is assumed to be infinite, and the 
QWR size in the x (y) direction is a x (a y ). Here the x, y 
and z directions correspond to the [100], [010] and [001] 
crystallographic directions, respectively. We assume that 
the conduction and valence bands are decoupled, and the 
effect of the split-off band is neglected because the spin- 
orbit split-off energy in bulk GaAs is much larger than 
the energy gap between the subbands caused by the con- 
finement considered here. Then, based on the four-band 
Luttingcr-Kohn model^ the explicit matrix form of the 
4x4 bulk- hole Hamiltonian in the basis of spin-3/2 pro- 
jection (J z ) eigenstates with quantum numbers +|, +~, 
— h and — | can be written as^ 
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FIG. 1: Schematic geometry of the QWR. 



In these equations, too denotes the free electron mass, 71, 
72 and 73 are the Luttinger parameters, £ is the electric 
field, and Ji are spin-3/2 angular momentum matrices. 4 - 
HEv8v i s the SOC arising from the structure inversion 
asymmetry (SIA) and H% v&v is the SOC from the bulk 
inversion asymmetry (BIA). These two terms turn out 
to be one or two orders of magnitude smaller than the 
intrinsic SOC from the four-band Luttinger-Kohn Hamil- 
tonian [the first term in Eq. ([!])] . This can be seen from 
Appendix A where we present a comparison between spin 
splittings due to the SIA and BIA and the splitting from 
the intrinsic SOC. Moreover, from the first term in Eq. 

, one can see that the LH spin-up states can be directly 
mixed with the HH states by S and R, but the mixing be- 
tween LH spin-up states and LH spin-down states has to 
be mediated by the HH states. All the mixing is related 
to the confinement. When the confinement decreases, the 
mixing increases due to the decrease of the energy gap 
between the LH and HH states. 

We construct the KSBEs by using the nonequilibrium 
Green function method as follows4£ 



Pk + pk\coh + pk\scatt = 0. 



(8) 



Here pk represents a single-particle density matrix of 
holes with wave vector k along the z-axis. One can 
project pk in the collinear spin space which is constructed 
by basis {s}, with {s} obtained from the eigenfunctions 
of the diagonal part of Hh(k). \s) = \m,n)\cr) with 
(r\m,n) = -?M= sin(^) sin (^) e lkz and \a) stand- 

' ' ^a x a y \ a y ' v a x > I ' 

ing for the eigenstates of J z . Then the matrix elements 
in the collinear spin space p\ S2 is written as p c k S2 = 
(si\pk\s2)- Here the superscript "c" denotes the quantum 
number distinguishing states in the collinear spin space. 
One can also project p k in the "helix" spin space which is 
constructed by basis {77} with 77 being the eigenfunctions 
ofF fc (fc): 



H h {k)\ri) =E nik \ri). 



(9) 



This basis function is a mixture of LH and HH states and 
is k dependent. Then the matrix elements in the helix 
spin space p%, n , n > can be written as p\ r) „, = (r}\pkW), 
with the superscript u h" denoting the helix spin space. 
The density matrix in the helix spin space can be trans- 
formed from that in the collinear one by a unitary trans- 
formation: p\ = Uj.p'f.Uk, where Uk{i,a) — rf a (k) with 
rf a {k) being the ith element of the ath eigenvector after 
the diagonalization of Hh(k). 
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In this paper we project the density matrix in the helix 
spin spaced and then the coherent terms can be written 

as 



Pk \coh 
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Iq is the form factor in the collinear spin space with wave 
vector Q = {q xi q yi q)- The first term in Eq. (fTU)) is the 
Coulomb Hartree-Fock term, and the second term is the 
contribution from the intrinsic SOC from the Luttingcr- 
Kohn Hamiltonian. Iq can be written as /q. Si .s 2 — 
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For small spin polarization, the contribution from the Hartree-Fock term in the coherent term is negligibl e) 1 and 
the spin precession is determined by the SOC, p ki „„/\coh = ^iPk^rf (^,k — E v i tk ), which is proportional to the energy 
gap between 77 and 77' subbands. 

The scattering terms include the hole-nonmagnetic-impurity, hole-phonon and hole-hole Coulomb scatterings. In 
the helix spin space, The scattering terms are given by 
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in which T k _ ri {i,j) = S n iS V j. Vq in Eq. (p~2|) reads 
Vq = 47re 2 /[Ko(c/ 2 + q 2 + k 2 )], with ko representing the 

static dielectric constant and k 2 = 4ire 2 Ni l /(kBTKo<i 2 ) 
standing for the Debye screening constant. N{ in Eq. (fTSj) 
is the impurity density and \Uq\ 2 = {AnZie 2 /[no(q 2 + 
q 2 + k 2 )]} 2 is the impurity potential with Zi standing 
for the charge number of the impurity. |Mq a| 2 and 
Nq^\ — [exp(cjQ .x/ksT) — 1] _1 are the matrix element 
of the hole-phonon interaction and the Bose distribution 
function with phonon energy spectrum ojq^x at phonon 
mode A and wave vector Q, respectively. Here the 
hole-phonon scattering includes the hole-LO-phonon and 
hole-AC-phonon scatterings with the explicit expressions 
of |A/q,a| 2 can be found in Refs.HIl 

It is noted that in the scattering terms Eq. (fT2")l . 
the energy spectra E Vtk are from Eq. © with full 



SOC included. As discussed by Cheng and Wu, this 
spectrum leads to the so called helix statistics in the 
equilibrium^ i.e., the Fermi distribution with SOC in- 
cluded in the energy spectrum. In most of our previ- 
ous works ) 35 i 36 i2L2£ i^ji-JiJ-L-iil the energy spectra in the 
scattering terms do not include the SOC, i.e., the en- 
ergy spectra from the Hamiltonian without the SOC. 
The corresponding equilibrium statistics is referred to as 
the collinear statistics^ It has been demonstrated that 
when the SOC is weak, the collinear statistics is good 
enough. However, the SOC for holes in QWR system 
can be very strong. Therefore, it is important to adopt 
the helix statistics in this investigation. 

Finally we comment on the reason we solve the KSBEs 
in the helix spin space p\. This is because the numerical 
calculation becomes faster in the helix spin space. This 
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can be understood from the fact that even the lowest 
helix subband is a mixture of many collinear basis states 
{s} due to the strong SOC. Therefore, a large number of 
basis states have to be included in the calculation if we 
project the KSBE in the collinear spin space. 

III. NUMERICAL RESULTS 

We first solve Eq. © to obtain the subband structure. 
In Fig. 2 we show six typical energy spectra for different 
confinements. Each subband is denoted as 1+ (I—) if the 
dominant spin component is the spin-up (-down) state. 
One can see from Fig. 2 that 1+ and 1— are very close to 
each other, so are the subbands 2±. The spin-splitting 
between them is mainly caused by the SOC arising from 
BIA, for that the spin-splitting caused by the SOC arising 
from SIA is three orders of magnitude smaller than the 
diagonal terms in Eq. JT]) and can not be seen in Fig. 2. 
The spin-splitting caused by the BIA is proportional to 
(P,~ — Py)-> which disappears when the confinement in 
x and y directions are symmetrical. Therefore, Z± are 
almost degenerate when a x = a y = 10 nm. If one ex- 
cludes the SOC from the BIA and SIA, i± are always 
degenerate because of the Kramers degeneracy. One also 
observes that when a x gets larger, the subbands are closer 
to each other. Especially, in the case of a x = a y = 10 nm, 
there are anticrossing points due to the HH-LH mixing in 
the Luttinger Hamiltonian. When a x keeps on increas- 
ing, the anticrossing point at small k between the 1± and 
2± gradually disappears. However, at large k region, the 
lowest two subbands become very close to each other. 
These will lead to significant effect on SRT. 

In order to show the situation of hole's population 
in these subbands clearly, We introduce a quantity (E), 
with 

= E, /_ ~ dk(pl +A+ ri il _ tl _)(E l+ , k + £ t _, fc ) 

2E I /_ + :*(pi, i+ , i+ -Pt-, I -) 

(13) 

to represent the energy region where spin precession and 
relaxation between the + and — bands mainly take place. 
Due to the small spin polarization, (E) is approximately 
equal to the Fermi energy at very low temperature. In 
Fig. 2 we plot (E) for N h = 4x 10 5 cm" 1 and 2x 10 6 cm" 1 
at T = 20 K. It is seen that (E) only intersects with the 
1± and 2± subbands, which means holes are mainly pop- 
ulated in 1± and 2± subbands^ Therefore, only the 1± 
and 2± subbands are taken into account in the present 
investigation. Higher subbands should be included if 
one considers higher hole density, higher temperature, or 
larger QWR size. It is further noticed that the dominant 
spin component in 1+ (1—) state is the spin- up (spin- 
down) LH state. At small k, the spin-up (spin-down) LH 
admixture remains at more than 90 %. Moreover, the 
HH-LH mixing in 2± subbands is much stronger. 

After the energy spectrum is obtained, we numerically 
solve the KSBEs and obtain the temporal evolution of the 
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FIG. 2: Typical energy spectra for (a) a x = 6 nm; (b) a x = 
8 nm; (c) a x = 10 nm; (d) a x = 12 nm; (e) a x = 15 nm; 
and (f) a x = 20 nm. a v = 10 nm. (E) at T = 20 K is also 
plotted: solid line for Nh = 4x 10 5 cm -1 and dashed line for 
N h = 2 x 10 6 cm" 1 . 

hole density matrix p^(t) in helix spin space. Then we 
project pfe(t) back into the collinear spin space p%(t) and 
obtain the temporal evolution of the spin polarization 

Sl{t)=Tv[pl{t)3], (14) 

in which J is the operator for spin-3/2 angu- 
lar momentum, written as a matrix in the ba- 
sis of z-projection eigenstates with eigenvalues m — 
+3/2, +1/2, -1/2, -3/2. We include the hole-phonon 
and hole-hole scatterings throughout our computation. 
The material parameters of GaAs in our calculation are 
the same as those used in Refs. 44.451 The initial con- 
dition at t = is set to be spin polarized with a small 
initial spin polarization P which is given, in the helix 
spin space, by P = (JVi+ + N 2 + - Nx- - N 2 -)/N h where 
Nh is the total hole density. Therefore, we have initial 
spin polarization along all directions in the collinear spin 
space. Then as discussed in the previous papers^ the 
SRT t can be defined by the slope of the envelope of the 
spin polarization along the z-axis: 

k 

A. Spin relaxation mechanisms 

There are three mechanisms leading to spin relaxation. 
First, the spin-flip scattering, which includes the scatter- 
ing between 1+ and I— subbands and the scattering be- 
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tween 1+ and I'— subbands (I ^ I'), can cause spin relax- 
ation. The SRT decreases with the spin-flip scattering, 
with the scattering strength being proportional to the 
spin mixing of the helix subbands. Second, because of the 
coherent term p\\ co h-, there is a spin precession between 
different subbands. The frequency of this spin precession 
depends on k and this dependence serves as inhomoge- 
neous broadening. As shown in Refs. 35.36.37. 3 S'l l 391140 . 
in the presence of the inhomogeneous broadening, even 
the spin-conserving scattering can cause irreversible spin 
relaxation. As a result, the spin-conserving scattering, 
i.e., the scattering between 1+ and l'+ and the scatter- 
ing between I— and I'—, can cause spin relaxation along 
with the inhomogeneous broadening. At last, the spin- 
flip scattering along with the inhomogeneous broadening 
can also cause an additional spin relaxation. 

It is seen from Fig. 2(a) that when Nh = 4 x 10 5 cm -1 
and a x = 6 nm, (E) only intersects with the 1± subbands 
and is far away from the 2± subbands. Therefore, holes 
populate the 1± subbands only. As pointed out before, 
the coherent term p\ 1+ 1 _ \ co h is proportional to — 
Ei^ k). As holes are only populating states in the small k 
region where the spin splitting between 1± is negligible, 
the spin precession between these two states, and thus the 
inhomogeneous broadening, is very small. Consequently 
the main spin-relaxation mechanism is due to the spin- 
flip scattering, i.e., the scattering between 1± subbands. 

In the case of larger a x and Nh as shown in Fig. 2(c)- 
(f), where (E) is close to or intersects with the 2± sub- 
bands, holes populate both the 1± and 2± subbands. 
The spin-flip scattering here includes the scattering be- 
tween 1± states, the scattering between 2± states and 
the spin-flip scattering between 1± and 2± subbands. 
This spin-flip scattering is still found to be the main spin 
relaxation mechanism. Besides, differing from the case 
of Fig. 2(a), the coherent term p\ 1± 2± | co /j is propor- 
tional to the energy gap between 1± and 2±, and it is 
much larger than p\ 1+ 1 _\ co h- As a result, there is a 
much stronger spin precession between 1± and 2± sub- 
bands with a frequency depending on fc, and the inho- 
mogeneous broadening caused by this precession along 
with both the spin-conserving scattering and the spin- 
flip scattering can make a considerable contribution to 
the spin relaxation. 



B. Wire width dependence of the SRT 

In Fig. 3 we plot the SRT as a function of the QWR 
width in x direction, a Xl for various temperatures. Here 
a y = 10 nm, the hole density is taken to be Nh = 
4 x 10 5 cm" 1 in Fig. 3(a) and N h = 2 X 10 6 cm" 1 in 
Fig. 3(b). It is seen from Fig. 3(a) that in the case 
of low density and T = 20 K, the SRT first decreases 
with a x when a x < 10 nm, then increases with a x when 
10 nm < a x < 14 nm, and finally decreases with a x when 
a x > 14 nm. To understand this behavior, let's look at 
the energy spectra for a wire with a x = 6 nm and 10 nm 
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FIG. 3: SRT r vs. the QWR width in x direction a x for 
(a) Nh = 4 x 10 5 cm -1 at different temperatures and (b) 
N h = 2x 10 6 cm" 1 at T = 20 K. a v = 10 nm. 



in Fig. 2. When the wire width increases, the energy gap 
between the 1± and 2± becomes smaller, and the spin 
mixing in the helix subbands increases. Therefore, the 
contribution from all of the spin-flip scattering increases, 
and the SRT decreases with increasing a x . From this 
point of view, one can expect a minimum of SRT when 
all of the spin mixing reaches a maximum at the Fermi 
surface as we can approximately make the assumption 
that the spin relaxation occurs mainly around the Fermi 
surface. In order to show this effect, let's look at (E) 
at T = 20 K for N h = 4 X 10 5 cm" 1 in Fig. 2(c). One 
can see that in the case of a x — a y = 10 nm, the lowest 
two subbands have an anticrossing and (E) is very close 
to the anticrossing point. From our calculation, we find 
that all of the spin mixing, including the HH-LH mixing, 
the mixing between the LH up states and the LH down 
states, and the mixing between the HH up states and the 
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HH down states, reaches a maximum because of the an- 
ticrossing. As a result, this anticrossing point leads to a 
strong spin-flip scattering and accounts for the minimum 
of SRT in Fig. 3(a) at T = 20 K as expected. When a x 
keeps on increasing, (E) will move into the larger-fc re- 
gion and the anticrossing point will disappear gradually 
as shown in Fig. 2, and the energy gap between the low- 
est two helix subbands at the Fermi surface gets larger. 
Therefore, the spin mixing at the Fermi surface becomes 
smaller, and the SRT will slightly increase with a x . How- 
ever, when a x > 14 nm, the effect of reducing the energy 
gap between the lowest two subbands and increasing the 
spin mixing are more important and the SRT decreases 
with a x again. 

We also plot the SRT as a function of a x at 100 K and 
300 K in Fig. 3(a). One finds that the SRT decreases 
with temperature. This could be understood as follows: 
firstly, when the temperature increases, holes are popu- 
lating the higher-fc states, for which all of the spin mix- 
ings in the 1± and 2± subbands are stronger; secondly, 
the strength of total scattering is enhanced. Both effects 
increase the contribution of spin-flip scattering and speed 
up the spin relaxation. 

It is seen from Fig. 2(f) that the lowest two subbands 
are very close to each other in large k region. There- 
fore, one can expect a very short SRT when the Fermi 
surface enters this region as the spin mixing here is very 
large and the contribution from the spin-flip scattering 
could be very strong. In order to show this effect, we 
take the hole density to be JVj = 2x 10 6 cm -1 to place 
the Fermi surface into the larger k region in Fig. 3(b). 
It is seen from the figure that the SRT decreases mono- 
tonically with a x . This could be easily understood from 
the fact that the LH-HH mixing increases with a x due to 
the decrease of the energy gap between the LH and HH 
states. In the case of a x — 20 nm, the SRT is nearly two 
orders of magnitude smaller than the SRT at a x = 6 nm 
as expected. 

In Fig. 4 we plot the SRT as a function of a x with dif- 
ferent a y at N h = 4 x 10 5 cm" 1 and T = 100 K. It is 
seen that, when a y = 5 nm, the SRT decreases monoton- 
ically with a x . This is because when the confinement is 
strong, there is no anticrossing point between the lowest 
two subbands in the region where holes are distributed. 
As a result, the SRT decreases with a x because of the 
effect of reducing the energy gap between the lowest two 
subbands. When a y is increased, the anticrossing point 
appears and the SRT shows a minimum similar to the 
case shown in Fig. 3(a). 



C. Hole density and temperature dependence of 
SRT 

Now we turn to study the hole-density dependence of 
the SRT at different temperatures and confinement sizes. 
In Fig. 5(a) we plot the SRT as a function of Nh at var- 
ious temperatures and a x = a y = 6 nm. From the figure 
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FIG. 4: SRT r vs. the QWR width in x direction at different 
a y at N h = 4 x 10 5 cm" 1 and T = 100 K. 



one can see that, when T < 50 K, the SRT first increases 
then decreases with A^. To understand this behavior, 
let's look at the energy spectrum for a x = a y = 6 nm 
shown in Fig. 5(c). The dashed line in Fig. 5(c) repre- 
sents (E) for N h = 12 x 10 5 cm" 1 at T = 20 K. One 
can see that the energy gap between 1± and 2± is large 
because of the small QWR size. When the tempera- 
ture is low, (E) only intersects with the 1± subbands 
and is far away from the 2± subbands. Therefore, holes 
populate the 1± subbands only. The main spin relax- 
ation mechanism is from the spin-flip scattering, i.e., the 
scattering between 1± subbands. Furthermore, it can 
be seen that the 2± subbands have a maximum at the 
wavenumber where the 1± subbands intersect with (E) 
for N h = 12 x 10 5 cm" 1 at T = 20 K. In the region 
where k is smaller than the intersection point of (E) and 
the 1± subbands, the energy gap between 1± and 2± 
subbands increases, and our calculation shows that the 
spin mixing of the 1± states decreases. Therefore, when 
T = 20 K and N h < 12 x 10 5 cm" 1 , the spin mixing 
at the Fermi surface decreases with increase of Nh, and 
the SRT increases with Nh because of the decrease of the 
spin-flip scattering. In the region where k is larger than 
the intersection point, the energy gap between 1± and 2± 
subbands decreases, and the spin mixing of the 1± states 
increases. As a result, the SRT increases with Nh when 
N h > 12 x 10 5 cm" 1 . The case of T = 50 K is similar to 
that of T — 20 K, but the holes are distributed in a wider 
k region and reach the maximum of the 2± subbands at 
a smaller Nh- Therefore, the SRT begins to decrease at 
N h = 8 x 10 5 cm" 1 . One can further see that the SRT 
for T = 20 K is smaller than that for T = 50 K at small 
Nh but larger than it at large Nh- This can be easily un- 
derstood as when the Fermi wavevector is smaller than 
the wavevector where the maximum of the 2± subbands 
occurs, the spin mixing decreases with T and the SRT 
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FIG. 5: SRT vs. the hole density at different temperatures, 
(a) a x — a y — 6 nm; (b) a x = a y — 10 nm. The energy 
spectrum for a x = a y — 6 nm is shown in (c). 



increases with T. Otherwise, the SRT decreases with T. 

We also plot the SRT as a function of Nh at 100 K and 
300 K in Fig. 5(a). One finds that the SRT decreases 
with Nh- This can be understood as follows: firstly, holes 
are populated at high-fc states that are larger than the 
wavevector where the maximum of the 2± subbands oc- 
curs, and the spin mixing increases with Nh', secondly, 
due to the larger temperature, the holes are also dis- 
tributed in 2± states. Therefore, the spin-flip scattering 
includes not only the scattering between 1± subbands, 
but also the inter-subband spin- flip scattering, i.e., the 
spin-flip scattering between 1± and 2± subbands. The 
strength of this inter-subband spin-flip scattering is en- 
hanced with the increase of Nh because of the increase 
of the hole population in 2± states. Both effects increase 
the contribution of spin-flip scattering and boost the spin 
relaxation. 

The results shown in Fig. 5(a) are quite different com- 
pared with those of LHs in quantum wells with small 
well width where the SRT decreases monotonically with 
Nh at low temperature but increases with Nh at high 
temperature^ The difference originates from the fact 
that the energy spectrum of the QWR is modulated dra- 
matically by the QWR size, and one can modulate the 
spin mixing strength by changing the region where holes 
are distributed. In order to better show this modula- 
tion, we also plot the results of a x — a y = 10 nm in Fig. 
5(b). In this situation the energy gap between the low- 
est two subbands gets smaller as shown in Fig. 2(c), and 
consequently both the spin mixing and the strength of 
the inter-subband spin-flip scattering increase with Nh 
like the case of T > 100 K in Fig. 5(a). Therefore, 
the SRT decreases with Nh as expected. The case of 
Nh > 5 x 10 5 cm -1 is not calculated in Fig. 5(b) as higher 
subbands should be included when Nh > 5 x 10 5 cm -1 
and T > 200 K, while only the 1± and 2± subbands are 
taken into account in our investigation. 

To see more detail of how the temperature affects the 
spin relaxation, we plot in Fig. 6 the SRT as a function 
of T with the QWR size of a x = a y = 6 nm. When 
Nh = 4 x 10 5 m _1 , the SRT first increases then decreases 
with T for the reason mentioned in the paragraph above. 
When Nh = 2 x 10 6 m _1 , one can see a fast decrease 
of the SRT around T = 100 K. To understand this be- 
havior, we plot the hole distribution in the 2± subbands 
in Fig. 6, from which one can see a fast increase of pop- 
ulation in the 2± subbands around T = 100 K, as the 
energy scale of the gap between the 1± and 2± subbands 
is close to fcflT of T = 100 K. This fast increase of the 
hole occupation in the 2± subbands leads to an increase 
of inter-subband spin-flip scattering which accounts for 
the fast decrease of the SRT around T = 100 K. To fur- 
ther reveal the contribution of inter-subband scattering, 
we plot the results which exclude the inter-subband holc- 
phonon scattering and inter-subband hole-hole scattering 
as dashed curves in Fig. 6. One finds that the fast de- 
crease of the SRT around T = 100 K disappears. We 
also plot the results without the coherent term but in- 
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FIG. 6: SRT r vs. temperature T. A: N h = 4 x 10 s m" 1 ; 
• : N h = 2 x 10 6 m" 1 . The QWR size is a x = a y = 6 nm. 
The solid curve is the result from the full calculation; the 
dashed curve is the result without the inter-subband hole- 
phonon scattering and inter-subband hole-hole scattering; the 
chain curve is the result without the coherent term. We also 
plot the hole distribution of the 2± subbands. Note the scale 
of the distribution is on the right side of the figure. 



eluding all the scattering as the chain curve. One can 
see that when T < 100 K and the holes are populated 
at the 1± subbands only, the chain curve coincides with 
the solid curve for that the coherent term only includes 
Pk i+ i-\coh which is negligible. When T > 100 K and 
holes populate both the 1± and 2± subbands, there is dif- 
ference between the chain curve and the solid curve which 
is due to the inhomogeneous broadening in the coherent 
term 1± 2 ±\coh- This inhomogeneous broadening, to- 
gether with the inter-subband scattering, can cause spin 
relaxation as discussed in subsection A. However, the dif- 
ference between the chain and the solid curves is small. 
This indicates that the contribution of this spin relax- 
ation mechanism is not as important as the contribution 
from the spin-flip scattering. 



D. Spin polarization dependence of SRT 

Finally, we investigate the initial spin polarization de- 
pendence of the spin relaxation. In Fig. 7 we plot the SRT 
as a function of hole density for both low and high spin 
polarizations with different QWR sizes at T = 100 K. It 
can be seen from the figure that the SRT of the case with 
high spin polarization is larger. This originates from the 
Hartree-Fock contribution of the hole-hole coulomb in- 
teraction, which serves as an effective magnetic field and 
can effectively reduce the spin relaxation at large spin 
polarizationi 38 i 42 




1 2 3 4 5 



N h (10 5 cm" 1 ) 

FIG. 7: SRT vs. the hole density N h at different QWR sizes 
and initial spin polarizations. The solid (dashed) curves are 
the results obtained for an initial spin polarization P = 2.5% 
(40%). ■: a x = a y = 6 nm; A: a x = a y = 10 nm. T = 100 K. 



IV. CONCLUSION 

In conclusion, we have investigated the spin relaxation 
of holes in p-type GaAs QWRs. The SRT is calculated 
by numerically solving the fully microscopic kinetic spin 
Bloch equations in the helix spin space. Differing from 
our previous works in n-type quantum-well and QWR 
systems2£i2&2L2&212id£ where the SOC is weak and the 
collinear statistics is good enough, the helix statistics is 
adopted in this investigation because of the strong SOC 
for holes in QWR system. Using this approach, we have 
studied in detail how the hole spin relaxation is affected 
by the wire size, the hole densities, temperature and the 
spin polarization. The confinement potential is assumed 
to be rectangular hard wall potential with infinite-depth 
throughout the paper. In real sample, the SRT can be 
quantitatively different from our results due to the dif- 
ferent confinements. However, the leading features such 
as the strong HH-LH mixing and the anticrossing points 
will be retained^ and the qualitative results will remain 
unchanged accordingly. 

We show that when holes are populated at the lowest 
helix subbands 1± only, the main spin-relaxation mech- 
anism is the spin-flip scattering which is proportional to 
the spin mixing of the helix subbands. When holes are 
populated in both 1± and 2± subbands, there are three 
mechanisms leading to spin relaxation: first, the bare 
spin-flip scattering; second, the spin-conserving scatter- 
ing along with the inhomogeneous broadening; and third, 
the spin-flip scattering along with the inhomogeneous 
broadening. However, the bare spin-flip scattering is still 
the dominant spin relaxation mechanism. 

The QWR size influences the SRT effectively because 
the spin mixing and the subband structure in QWRs de- 
pend strongly on the confinement. When the wire width 
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gets larger, the subbands are closer to each other and 
the spin mixing of the subbands gets larger, therefore, 
the SRT decreases. Especially, in the case of a x — a y = 
10 nm, there is an anticrossing point. If the Fermi surface 
happens to be close to this point, this anticrossing point 
leads to a strong spin-flip scattering which accounts for 
a minimum of the SRT. When a x = 20 nm, a y = 10 nm, 
the anticrossing point at small k between the 1± and 
2± disappears. However, at large fc, the lowest two sub- 
bands become very close to each other. If we take the 
hole density to be Nh — 2 x 10 6 cm -1 to place the Fermi 
surface at this large k region, the SRT is nearly two or- 
ders of magnitude smaller than the SRT at small wire 
width because the spin mixing here is very large and the 
contribution from the spin-flip scattering is very strong. 

The hole density influences the SRT by modulating the 
strength of spin mixing and the strength of inter-subband 
spin-flip scattering. In most of the cases we considered, 
the strength of spin mixing and the inter-subband scat- 
tering increase with Nh as the holes are populated in high 
k states. As a result, the SRT decreases with Nh- How- 
ever, when the confinement is very strong and the energy 
gap between 1± and 2± is large, there is a small region 
where the spin mixing decreases with k. If we choose a 
small Nh and a low temperature to make the holes be 
distributed in this small region only, one finds that the 
SRT increases with Nh because of the decreasing spin 
mixing. 

The influence of temperature on the SRT is similar to 
the case of the hole density dependence. The strength 
of both the spin mixing and the spin-flip scattering in- 
creases with T, and the SRT decreases with T. Espe- 
cially, when the energy scale of the gap between the 1± 
and 2± subbands is close to fc^T, there is a fast increase 
of the distribution on the 2± subbands, which leads to an 
increase of inter-subband spin-flip scattering and leads to 
a fast decrease of the SRT. This decrease of the SRT also 
proves that the inter-subband spin-flip scattering makes 
a marked contribution to the spin relaxation. We further 
show that the Hartree-Fock term increases with spin po- 
larization and can reduce the spin relaxation. 
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APPENDIX A: A COMPARISON OF BIA, SIA 
AND THE INTRINSIC SOC 

The SOC contributions for holes arising from the BIA 
and SIA can be obtained by quasi-degenerate pertur- 



bation theory (Lowdin partitioning) from the extended 
Kane model.— For an external electric field along x di- 
rection, the dominant terms of the SIA contribution can 
be written, in an explicit matrix notation, as follows^ 9 - 



n SvSv — 



SvSvc 
741 c x 



Ip 

2 V 



x/3 



iP y 



\ 



-ip 

2 r V 
lP z 







-ip z 

ip 

4^ 



v2 

2 

2 r v 



iPz 



J 

(Al) 

The coefficient -y%f v for GaAs is -1.462 x 1(T 19 e-m 2 
and is two or three orders of magnitude smaller than the 
off-diagonal terms in Eq. ([I]). 49 Besides, this term couples 
the two LH states directly while in Eq. JTj the two LH 
states can only mix with each other mediated by the HH 
states. However, this direct coupling is still very small 
compared to the intrinsic mixing due to the first term in 
Eq. ((T|). To show this, we study a simplified case includ- 
ing only the lowest eight collinear subbands of |1,1, a) 



and II, 2, a) with a 



±| and a 



±i and we use the 



second-order Lowdin partitioning to convert this 8x8 
matrix expanded by Eq. ([TJ to a block-diagonal form in 
which the off-diagonal matrix elements between |1, 1, ±i) 
and the other states are zero. Then the effective coupling 
between these two lowest LH states can be written as: 



H (2) i = (l,l,i|^|L2,|)(l,2,||fl|l,l,-i) 
, (l,l,||i2|l I 2,-f)(l,2, 



%1,A - ^1,2,- 

h 2 6472 73 k 
2m a(37i - 137 2 ) ' 



(A2) 



in which we assume a x = a y = a for simplicity. Then 
we compare this term to the coupling between |1, 1, ±5) 

8-»- - - 

contributed by SIA, and find — 



= 4.1 x 10 3 when 



a = 6 nm and £ x = 100 kV/cm (in experiments with 
quantum wells, the values of £ x are typically of the order 
of several kV/cm)^ Therefore, the contribution from 
SIA is very small. 

The SOC contribution from the BIA is also very 
small. The BIA coefficient b\f v in Eq. © is -8.193 x 
10 _29 eV-m 3 for GaAs^ 9 - and is one order of magnitude 
smaller than the intrinsic mixing due to the first term in 
Eq. (JTJ) . Furthermore, when only the lowest four collinear 
states (two for HHs and two for LHs) are included, one 
finds from Eq. ([7]) that only the third term is nonzero. 
However, the third term is diagonal and does not con- 
tribute any coupling. Therefore, the spin coupling due 
to the BIA contributions must be mediated by higher 
collinear states. 



10 



* Author to whom correspondence should be addressed; 
Electronic address: mwwu@ustc.edu .cn. I 

1 D. D. Awschalom, D. Loss and N. Samarth, Semiconductor 
Spintronics and Quantum Computation (Springer, Berlin, 
2002); I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. 
Phys. 76, 323 (2004); J. Fabian, A. Matos-Abiague, C. 
Ertler, P. Stano, and I. Zutic, acta physica slovaca 57, 565 
(2007), and references therein. 

2 S. Datta and B. Das, Appl. Phys. Lett. 56, 665 (1990). 

3 J. Schliemann, J. C. Egues, and D. Loss, Phys. Rev. Lett. 
90, 146801 (2003). 

4 J. Zhou, Q. W. Shi, and M. W. Wu, Appl. Phys. Lett. 
84, 365 (2004); K. Shen and M. W. Wu, Phys. Rev. B 77, 
193305 (2008). 

5 L. N. Pfeiffer, R. de Picciotto, K. W. West, K. W. Baldwin, 
and C. H. L. Quay, Appl. Phys. Lett. 87, 073111 (2005). 

6 R. Danneau, W. R. Clarke, O. Klochan, A. P. Micolich, 
A. R. Hamilton, M. Y. Simmons, M. Pepper, and D. A. 
Ritchie, Appl. Phys. Lett. 88, 012107 (2006). 

7 O. Klochan, W. R. Clarke, R. Danneau, A. P. Micolich, 
L. H. Ho, A. R. Hamilton, K. Muraki, and Y. Hirayama, 
Appl. Phys. Lett. 89, 092105 (2006). 

8 Y. A. Bychkov and E. Rashba, Sov. Phys. JETP Lett. 39, 
78 (1984). 

9 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

10 W. Hausler, Phys. Rev. B 63, 121310 (2001). 

11 S. Pramanik, S. Bandyopadhyay, and M. Cahay, Phys. 
Rev. B 76, 155325 (2007). 

12 S. Zhang, R. Liang, E. Zhang, L. Zhang, and Y. Liu, Phys. 
Rev. B 73, 155316 (2006). 

13 D. Csontos and U. Ziilicke, Phys. Rev. B 76, 073313 
(2007); Appl. Phys. Lett 92, 023108 (2008). 

14 R. Danneau, O. Klochan, W. R. Clarke, L. H. Ho, A. P. 
Micolich, M. Y. Simmons, A. R. Hamilton, M. Pepper, D. 
A. Ritchie, and U. Ziilicke, Phys. Rev. Lett. 97, 026403 
(2006). 

15 Y. Arakawa, T. Yamauchi, and J. N. Schulman, Phys. Rev. 
B 43, 4732 (1991). 

16 U. Bockelmann and G. Bastard, Phys. Rev. B 45, 1688 
(1992). 

17 F. Vouilloz, D. Y. Oberli, M.-A. Dupertuis, A. Gustafsson, 
F. Reinhardt, and E. Kapon, Phys. Rev. Lett. 78, 1580 
(1997). 

18 D. S. Citrin and Y.-C. Chang, Phys. Rev. B 40, 5507 
(1989). 

19 O. Stier and D. Bimberg, Phys. Rev. B 55, 7726 (1997). 

20 N. Shtinkov, P. Desjardins, R. A. Masut, and S. J. Vlaev, 
Phys. Rev. B 70, 155302 (2004). 

21 G. Goldoni and A. Fasolino, Phys. Rev. B 52, 14118 

(1995) ; G. Goldoni, F. Rossi, E. Molinari, A. Fasolino, 
R. Rinaldi, and R. Cingolani, Appl. Phys. Lett 69, 2965 

(1996) . 

22 E. I. Rashba and E. Ya. Sherman, Phys. Lett. A 129, 175 
(1988). 

23 O. Bleibaum, Phys. Rev. B 71, 235318 (2005). 

24 P. I. Tamborenea, M. A. Kuroda, and F. L. Bottesi, Phys. 
Rev. B 68, 245205 (2003). 

25 J. Kainz, U. Rossler, and R. Winkler, Phys. Rev. B 68, 



075322 (2003). 

26 W. H. Lau, J. T. Olesberg, and M. E. Flatte, Phys. Rev. 
B 64, 161301 (2001). 

27 F. X. Bronold, A. Saxena, and D. L. Smith, Phys. Rev. B 
70, 245210 (2004). 

28 A. A. Burkov and L. Balents, Phys. Rev. B 69, 245312 

(2004) . 

29 S. Saikin, M. Shen, M.-C. Cheng, and V. Privman, J. Appl. 
Phys. 94, 1769 (2003). 

30 M. Ohno and K. Yoh, Phys. Rev. B 75, 241308 (2007). 

31 S. Pramanik, S. Bandyopadhyay, and M. Cahay, Phys. 
Rev. B 68, 075313 (2003); ibid. 73, 125309 (2006); P. 
Upadhyaya, S. Pramanik, S. Bandyopadhyay, and M. Ca- 
hay, ibid. 77, 045306 (2008). 

32 A. Dyson and B. K. Ridley, Phys. Rev. B 72, 045326 

(2005) . 

33 P. Schwab, M. Dzierzawa, C. Gorini, and R. Raimondi, 
Phys. Rev. B 74, 155316 (2006). 

34 T. Sogawa, H. Ando, and S. Ando, Phy. Rev. B 58, 15652 
(1998). 

35 M. W. Wu and H. Metiu, Phys. Rev. B 61, 2945 (2000). 

36 M. W. Wu and C. Z. Ning, Eur. Phys. J. B 18, 373 (2000). 

37 M. W. Wu, J. Phys. Soc. Jpn. 70, 2195 (2001). 

38 M. Q. Weng and M. W. Wu, Phys. Rev. B 68, 075312 
(2003); ibid. 70, 195318 (2004). 

39 M. Q. Weng, M. W. Wu, and L. Jiang, Phys. Rev. B 69, 
245320 (2004). 

40 For brief review, see M. W. Wu, M. Q. Weng, and J. L. 
Cheng, in Physics, Chemistry and Application of Nanos- 
tructures: Reviews and Short Notes to Nanomeeting 2007, 
eds. V. E. Borisenko, V. S. Gurin, and S. V. Gaponenko 
(World Scientific, Singapore, 2007) pp. 14, and references 
therein. 

41 J. L. Cheng and M. W. Wu, J. Appl. Phys. 99, 083704 

(2006) . 

42 D. Stich, J. Zhou, T. Korn, R. Schulz, D. Schuh, W. 
Wegscheider, M. W. Wu, and C. Schiiller, Phys. Rev. Lett. 
98, 176401 (2007); Phys. Rev. B 76, 205310 ' (2007). 

43 J. Zhou and M. W. Wu, Phys. Rev. B 77, 075318 (2008). 

44 J. Zhou, J. L. Cheng, and M. W. Wu, Phys. Rev. B 75, 
045305 (2007). 

45 C. Lii, J. L. Cheng, and M. W. Wu, Phys. Rev. B 73, 
125314 (2006). 

46 J. L. Cheng, M. Q. Weng, and M. W. Wu, Solid State 
Commun. 128, 365 (2003). 

47 T. Dietl, H. Ohno, F. Matsukura, J. Cibert, and D. Fer- 
rand, Science 287, 1019 (2000). 

48 J. Luttinger, Phys. Rev. 102, 1030 (1956). 

49 R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, Berlin, 
2003). 

50 From our calculation, we find that in the case of a x = 
a y = 20 nm, N h = 2 x 10 6 cm" 1 and T = 20 K, (E) is 
close to 3± states and the total population in 3± states is 
approximately 8 %. 

51 R. Winkler, S. J. Papadakis, E. P. De Poortere, and M. 
Shayegan, Phys. Rev. Lett. 84, 713 (2000). 



